Calculation of microscopic nuclear level densities based on covariant density functional theory 
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A microscopic method for calculating nuclear level density (NLD) is developed, based on the covariant den- 
sity functional theory (CDFT). The particle-hole state density is calculated by combinatorial method using the 
single-particle levels schemes obtained from the covariant density functional theory, and then the level densities 
are obtained by taking into account collective effects such as vibration and rotation. Our results are compared 
with those from other NLD models, including phenomenological, microstatistical as well as non-relativistic 
HFB combinatorial models. The comparison suggests that besides some small deviations from different NLD 
models, the general trends among these models are basically the same. In addition, the NLDs of the CDFT com- 
binatorial method with normalization are compared with experimental data, including the observed cumulative 
number of levels at low excitation energy and the measured NLDs. Compared with the existing experimental 
data, the CDFT combinatorial method can give reasonable results. 
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I. INTRODUCTION 


Nuclear level density (NLD) is the basic physics input for 
nuclear reactions. It is the key ingredient for the calculation 
of reaction cross sections relevant for nucleosynthesis [1-3]. 
The study of NLDs can date back to 1930’s with Bethe’s pio- 
neering work [4]. Since then, many theoretical models, such 
as the Back-Shifted Fermi gas model (BFM) [5], the Compos- 
ite Gilbert-Cameron model (CTM) [6] and the Generalised 
superfluid model (GSM) [7], have been adopted for the NLD 
study. These phenomenological models are widely used in 
nuclear reaction calculations. To adjust the parameters, the 
phenomenological models rely more or less on experimental 
data, however, the experimental data is limited, especially for 
nuclei far from the 6-stability line [8]. To deal with such dif- 
ficulties, the microscopic method has then been developed. 

Over the last decades, various microscopic models for 
NLD have been proposed, from the equidistant spacing 
model [9-12], the shell-model Monte Carlo method [13- 
17], the spectral distribution calculation [18-20], to the in- 
dependent particle model at finite temperature [21-24] and 
the micro-statistical methods [25-28]. Especially in the last 
two decades, microscopic method based on the Hartree-Fock- 
Bogoliubov (HFB) plus combinatorial model [8] has been 
well developed. The combinatorial method can compete with 
statistical methods in reproducing experimental data, and pro- 
vide the energy-, spin- and parity-dependent NLDs which 
are beyond the reach of statistical methods [29]. The non- 
relativistic Hartree-Fock-Bogoliubov combinatorial methods 
based on the Skyrme and Gogny effective interactions have 
successfully reproduced the NLDs for various nuclei [29, 30] 
and been applied to various astrophysical reactions. The ac- 
curacy of NLD is related to the basic information of nu- 
clear structure, such as single-particle levels, deformation and 
binding energy. In recent years, the covariant (relativistic) 
density functional theory has attracted considerable attention 
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in the nuclear physics community on account of its success- 
ful description of the complex nuclear structure and reaction 
dynamics [31-36]. For instance, it can reproduce well the 
isotopic shifts in the Pb isotopes [37] and naturally give the 
origin of the pseudospin and spin symmetries in the antinu- 
cleon spectrum [38], as well as provide a good description of 
the nuclear magnetic moments [39, 40]. Recently, a micro- 
statistical method based on CDFT has been developed to de- 
scribe NLDs [41]. The model is applied to the calculation of 
NLDs of 94:96-98Mo, 106,108pd, etc, and the NLDs are in very 
good agreement with experimental data over the entire energy 
range of measured values [41]. While the microstatistical 
method can calculate only energy-dependent NLDs, the com- 
binatorial method can calculate the energy-, spin- and parity- 
dependent NLDs. Therefore, it is a meaningful attempt to 
calculate the NLDs using CDFT plus combinatorial method. 

The theoretical framework and methods are introduced in 
Sec II. The NLDs calculated based on CDFT combinatorial 
method are compared with other NLD predictions and exper- 
imental data in Sec III. Conclusions and prospects are finally 
drawn in Sec IV. 


I. THEORETICAL FRAMEWORK 


Covariant density functional theory starts from a La- 
grangian, and the corresponding Kohn-Sham equations have 
the form of a Dirac equation with effective fields S and V de- 
rived from this Lagrangian [36, 40, 42-44]. Specifically, the 
nucleons in the nucleus are described as Dirac particles mov- 
ing in the average potential field given by the meson and pho- 
ton fields, interacting with each other through the exchange 
of meson and photon. By solving the Dirac equation: 


(1) 


where e; is the single-particle energy, and that’s what we need 
to calculate the NLDs. The S and V are the relativistic scalar 
field, and the time-like component of vector field. 

After obtaining the energy €, spin projection m and parity p 
of the single-particle levels through CDFT, the level informa- 
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tion is substituted into the generating function defined in the 
combinatorial method [8] in order to obtain the particle-hole 
state density p;, and the generating function Z reads 

4 Ir 


_ kak 
2 (a1, £2, £3, %4) = II [[a + pe yei t™ ). 
k=1i=1 


(2) 


This generating function is a straightforward generalization 
of that used previously in Ref. [12] to account not only for 
energy £" but also for spin projection m* and parity p*. Z 
can be expanded in powers of zy writing 


4 
2 (m1, 22, 03,04) = X Fuly, t) |] a’, (3) 
N k=1 


the symbol N denoting again any integers combination 
(Ni, N2, N3, N4). The function F(y,t) can be expanded 
into powers of y and t writing 


NGao=>) > 5 Cyw (U, M, P)y”t™, 


U M P=-1,+1 


(4) 


where the U, M and P are the excitation energy, spin projec- 
tion and parity. The coefficients Cy (U, M, P) are the num- 
bers of solutions. Through the coefficients Cy (U, M, P), the 
simplest definition for particle-hole state density p; reads 


= Sect M,P), 
E0 

but the state density p; turn out to be strongly dependent on 
any unit energy &o. Therefore, another method suggested by 
Williams [45] is employed in this paper to limit the discretiza- 
tion effects as explained below. 

Summing all the Cw up to a given excitation energy U, 
we first obtain the cumulated number of states Ny (U, M, P) 
which represents the number of particle-hole states with ex- 
citation energy E such that 0 < E < U. The particle-hole 
state density defined as 


dNy (U, M, P) 
p= dU 

The particle-hole state density p; calculated is only related 

to the particle-hole excitation. Two special collective effects 

need to be taken into account in order to obtain the level den- 

sity, including rotational, and vibrational enhancement. If the 

nucleus under consideration displays spherical symmetry, the 

intrinsic and laboratory frames coincide, and the level density 
is trivially obtained through the relation [46] 


Pspn (U, M, P) = 9(U,M = J, P) — pi(U,M = J + 1, P). 

(7) 

For deformed nuclei, within the axial symmetry hypothe- 
sis, the NLD after accounting for rotational effects reads 


. (6) 
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K=-J,KZ0 
+ (6(Jeven)O(p=+) pi(U = EZE, 0, P) 
+ 5(Jodd) 9(p=—) Pi (U = Ee 0, P)). 
(8) 
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Factor 1/2 accounts for the fact that in mirror axially sym- 
metric nuclei, the intrinsic states with spin projections +% 
or —K give rise to the same rotational levels. Moreover, in 
the second term of the summation, the symbol 6(,,.) (defined 
by O(a) = 1 if x holds true, and 0 otherwise) restricts the 
rotational bands built on intrinsic states with spin projection 
K = 0 and parity P to the levels sequences 0, 2, 4,... for 
P = + and 1, 3,5,... for P = —. Finally, the rotational 
energy is obtained with the well-known expression [47] 


gIK _ J(J+1)-K? 
rot `~ 27 ’ 
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where 7, is the moment of inertia of a nucleus rotating 
around an axis perpendicular to the symmetry axis. In this ar- 


(9) 


ticle 7, is approximated by the rigid-body value 77'*"" rigid 
which reads 

= 2 5 

[eit = ZmR?(1+ 4/6" 10 


for an ellipsoidal shape with axial quadrupole deformation 
parameter 6. The vibration enhancement is approximated by 
equation [48] 

Kyi = exp[6S = ($SU/T)], (11) 
where S is the entropy, U is the excitation energy, and T is 
the nuclear temperature. 

Finally, a phenomenological damping function is intro- 
duced to avoid sharp transitions between spherical and de- 
formed level densities affecting the NLD predictions. The 
expression of the NLD after introducing the damping func- 
tion can be written as [28] 


p(U, M, P) m [1 u faam(U)| vin Pspn(U, M, P) 


(12) 
7 faam(U)KvibPaet(U, M, P). 
The damping function faam is expressed as [30] 
1 
faam = (13) 


1 + exp (3 — 0.18)/0.038’ 


where 8 is the quadrupole deformation parameter. Param- 
eters 0.18 and 0.038 have been adjusted according to the ex- 
perimental data at the neutron separation energy Sn [49]. The 
expression of faam has been evolved in Ref. [30], which has 
been simplified to depend only on the nuclear deformation 
thus reducing the number of phenomenological parameters. 
This damping function is used to suppress discontinuities that 
occur between spherical and deformed NLDs and ensure a 
smooth shape change from deformed to spherical. 


Ill. RESULTS AND DISCUSSION 


In this section, we present our results of the NLDs from 
the combinatorial method based on CDFT and compare them 
with the results from other NLD models [50] and experi- 
ments [51, 52]. The effective meson-exchange interaction 


154 parameter PK1 [53] is adopted throughout the CDFT calcu- 
155 lations. For spherical CDFT calculations, we fix the box size 
156 Rpox = 20 fm, and step size Ar = 0.1 fm. In the present de- 
157 formed CDFT calculations, both the Dirac equation for nucle- 
158 ons and the Klein—Gordon equations for mesons are solved in 
159 an isotropic harmonic oscillator basis and a basis of 18 major 
160 oscillator shells is adopted. 


FT 
108 F112Cd (a) #62Dy (b) 
Alt + 
>= F = E 

oO 
S 10 F = EXP 4 
zx L — PKI 4 
10? L L ——NI3 | 
i + — DD-ME2 
10° ot PP PONY 
0 2 4 6 8 10 2 4 6 8 10 12 


E(MeV) E(MeV) 


Fig. 1. (Color online) Comparison of NLDs calculated using differ- 
ent effective interactions under CDFT with experimental data [51, 
52]. 
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As the combinatorial results rely on the single-particle lev- 
els properties. To understand the effect of using different 
effective interactions on the NLDs, the NLDs of !!?Cd and 
162Dy obtained with different CDFT effective interactions 
PK1 [53], NL3 [54], DD-ME2 [55] and DD-PC1 [56] are pre- 
sented for the sake of comparison in Fig. 1. For nuclei with 
small deformation, !!2Cd, the NLDs for four effective inter- 
actions are close for excitation energy above 4 MeV. How- 
ever, a significant difference is observed for excitation energy 
below 4 MeV. The difference comes from the fact that the 
single-particle energies near the vicinity of the Fermi level 
are very sensitive to the choice of effective interactions es- 
pecially when the Fermi level is near the proton or neutron 
shells. Meanwhile, CDFT calculations with chosen interac- 
tions do not reproduce the single-particle levels of the magic 
nuclei 13?Sn [57] very well, which explains the deviations 
of our results from the measurements in Fig.l. For the well 
deformed !®*Dy, the NLDs calculated with four effective in- 
teractions deviate from each other for the whole region of ex- 
citation energy, although the overall trend is consistent. This 
deviation at low excitation energy is mainly caused by differ- 
ences in the single-particle energies around the Fermi level. 
Meanwhile, the entropy S obtained from the four effective 
interactions is significantly different for the nuclei with large 
deformation, and ultimately leading to the deviation of NLDs 
at high excitation energy. Compared to the experimental data 
for !!2Cd [51] and 16?Dy [52], the deviation between the 
NLDs obtained by using PK1 effective interaction and the re- 
sults from other effective interactions is within a reasonable 
range. 

13 As the non-relativistic HFB combinatorial methods, in- 
194 Cluding Skyrme and Gogny interactions, have been widely 
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Fig. 2. (Color online) Comparison of NLDs calculated based on 
CDFT combinatorial method with results from HFB combinatorial 
methods and microstatistical method (HFB+Sta) [28-30]. The ex- 
perimental data are taken from Refs. [58-61]. 


195 used in NLD predictions [51, 62], the NLDs of spherical 
nuclei ( °°Co, !°7Sn, ?°8Pb ) and deformed nuclei ( °4Mo, 
95Mo, 166Er ) calculated from the present CDFT combinato- 
rial method are given in Fig. 2 and compared with them and 
the corresponding experimental data. 

In °°Co, the CDFT combinatorial method gives a lower 
NLD compared to the experimental data, while the Skyrme 
combinatorial method gives a higher NLD. Although the re- 
sult of Gogny combinatorial method is closer to the exper- 
imental data to a certain extent, there is a strong oscilla- 
tion. For !??Sn and 2°8Pb, the results of CDFT combinatorial 
method are similar to that of Skyrme combinatorial method, 
but the results of Gogny combinatorial method are signifi- 
cantly lower than the two methods. Compared with the ex- 
perimental data, the CDFT combinatorial method can better 
describe the NLD of ?0®Pb. The NLDs of *°Co, !??Sn and 
208Pb from the three combinatorial methods are more or less 
oscillatory, which is because the three nuclei are spherical nu- 
clei with highly degenerate single-particle levels, and there is 
a strong shell effect. 

For °*Mo, both the results of the CDFT and Skyrme com- 
binatorial methods are closer to the experimental data than 
the Gogny combinatorial method. For 9°Mo and !66Er, the 
NLDs obtained by the three combinatorial methods are simi- 
lar and slightly larger than the experimental data. It should 
be emphasized that the deformed nuclei °4*Mo, °°Mo and 
166Fr break single-particle level degeneracy, thus obtain a 
smoother NLDs as shown in Fig. 2. It is worth noting that 
the CDFT combinatorial method gives higher NLDs for °*Mo 
and °°Mo. One possible reason is that the CDFT combinato- 
rial method uses rigid-body value when considering rotational 
effects, which is not appropriate for soft nuclei [63]. 

In addition, the NLDs of °°Co and 9*Mo are well described 
228 by the microstatistical model [28] relative to the CDFT com- 
229 binatorial method, while the description of the NLDs of Mo 
230 and '66Er are similar between these two models. In compari- 
231 SON, microstatistical model gives relatively smooth results for 
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spherical and deformed nuclei, which do not well reflect shell 
effects. Therefore, it can be obtained that the CDFT combi- 
natorial method has a similar ability to describe NLDs as the 
other two combinatorial methods and microstatistical model. 
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Fig. 3. (Color online) Comparison of NLDs obtained from CDFT 
combinatorial method with other NLD predictions [50] and exper- 
imental data [51, 52, 59, 64]. The full asterisk corresponds to the 
experimental data at the neutron separation energy Sn [49]. 


At the same time, it’s also necessary to compare with phe- 
nomenological models (BFM, CTM and GSM) that have been 
widely applied. In particular, the phenomenological models 
can well describe the experimental data at the neutron sepa- 
ration by fitting experimental data. The NLDs of even-even 
nuclei (11?Cd, !6?Dy), odd-A nuclei (°'V, °’Mo, !!°Sn) and 
odd-odd nucleus (©°Co) calculated from the present CDFT 
combinatorial method are given in Fig. 3 and compared with 
results of phenomenological models and the experimental 
data. It can be seen from Fig. 3 that the CDFT combinatorial 
method has similar ability to the phenomenological models 
in reproducing the experimental data on neutron separation 
energy, especially the CDFT combinatorial method can de- 
scribe the NLD of 51V well. In addition, the overall descrip- 
tion of experimental data below neutron separation energy by 
the CDFT combinatorial method is also similar to that of the 
phenomenological methods. In conclusion, while there are 
some differences in the results for each NLD model, the over- 
all ability to describe experimental data is similar. 

The most extensive and reliable source of experimental 
information on NLD remains the s-wave neutron resonance 
spacings Do [49] and the observed low-energy excited lev- 
els [49]. To measure the dispersion between theoretical and 
experimental Do, the frms factor is defined as 


= pi,” 
apes i 
a — Sin? th 14 
f a| mop) > 0 


where Din(Dexp.) is the theoretical (experimental) resonance 
spacing and N. is the number of nuclei in the compila- 
tion. The results for the ratios of Di, and Dexp. are shown 
in Fig. 4, and the fims = 3.62, a total of 66 nuclei are 
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Fig. 4. The ration of CDFT combinatorial method (Dn) to the ex- 
perimental (Dexp.) s-wave neutron resonance spacings compiled in 
Ref. [49]. 


counted. Under the same calculation conditions, the fims of 
Gogny (D1S) interaction based on HFB plus combinatorial 
method is 7.25 [8]. At present, the result of CDFT combi- 
natorial method is bigger than frms = 1.8 deviation of the 
phenomenological back-shifted Fermi gas model [7]and the 
frms = 2.1 value obtained with microstatistical method [28]. 
But the HFB combinatorial based on Skyrme effective inter- 
action model improved pairing correlation [48] and collective 
effects, then get the firms = 2.3 [29]. In subsequent work, it is 
expected to improve the CDFT combinatorial method in the 
same way. 
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Fig. 5. (Color online) The NLDs obtained from CDFT combinato- 
rial method with (red) and without (black) normalization. The ex- 
perimental data are taken from Refs. [49, 59, 64]. 


When phenomenological NLDs are used in nuclear physics 
applications, such as nuclear data evaluation or accurate and 
reliable estimation of reaction cross sections, one often plays 
with the few parameters on which phenomenological expres- 
sions depend. The results of the combinatorial method also 
can normalize both the experimental level scheme at low ex- 


citation energy and the neutron resonance spacings at U = Sn 
in a way similar to what is usually done with analytical for- 
mulas. More precisely, the normalized level density can be 
obtained by expression [29] 


plU, J, P) = eV"? x p(U-6,5,P), 5) 
where the energy shift 6 is essentially extracted from the anal- 
ysis of the cumulative number of levels, and the a can be ob- 
tained through the expression 

293 Pth( Sn) x coven = Pexp. (Sn). (16) 
With such a normalization, the experimental low-lying states 
and the Do values can be reproduced reasonably well. 
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Fig. 6. Normalized parameters a (left) and 6 (right) values are plot- 
ted as a function of the atomic mass number. 


As an illustration, the variation of NLDs obtained by the 
CDFT combinatorial method before and after normalization 
is shown in Fig. 5. When the normalization is applied, the 
NLDs will pass through the experimental data at the neutron 
separation energy Sn. As shown in Fig. 5(a), if the theoret- 
ical calculation value is close to the experimental data at the 
neutron separation energy Sn, the results before and after nor- 
malization will not change much. In Fig. 5(b), the results will 
304 better agreement experimental data [64] after normalized, es- 
pecially above 2 MeV. This is because the coefficient e2 VU 3 
is related to the excitation energy U, and when U is small, the 
normalization has little effect. Some nuclear a and 6 values 
are shown in Fig. 6. 

Finally, the NLDs calculated based on the CDFT combina- 
torial method are compared with the observed low-energy ex- 
cited levels, which are the most extensive and reliable source 
of experimental information on NLDs [49]. The cumulative 
number of nuclear levels N (U) indicates the sum of the num- 
314 ber of all levels below some excitation energy U (including 
ais U) 
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where N(U, M, P) is the cumulative number of levels under 
spin M and parity P to the excitation energy E. The pre- 
dicted cumulative number of levels N(U) are compared with 
the experimental data [49] in Fig. 7, including light as well 
as heavy, spherical and deformed nuclei. Globally, the results 
of the CDFT combinatorial method are in reasonably good 
agreement with experimental data at lower excitation energy, 
especially for light nuclei. At high excitation energy, the the- 
oretically calculated cumulative number of nuclear levels is 
higher than the experimental value. NLDs increase rapidly 
with the increase of excitation energy, and the N(U) should 
also gradually increase at high excitation energy, but the ex- 
perimental levels of the cumulative number gradually stabi- 
lized, it may be limited by the experimental conditions and 
cannot completely count the number of levels. 

In Fig. 8, the CDFT combinatorial method predictions af- 
ter normalization are compared with the experimental data 
extracted by the Oslo group [51, 52, 59-61, 64-66]Jand the 
particle evaporation spectrum [67, 68]. The Oslo method is 
model-dependent. In order to extract the absolute value of 
the total level density from the measured data, the so-called 
experimental NLDs need to be normalized by the total level 
density at the neutron binding energy, which in turn is de- 
rived from the neutron resonance spacing. For a meaningful 
comparison between the CDFT combinatorial predictions and 
the Oslo group data, it is therefore important to normalize the 
NLDs of the CDFT combinatorial method to the level density 
value at U = $, considered by the Oslo group. As shown 
in Fig. 8, the results of the CDFT combinatorial method af- 
ter normalization agree well with the experimental data be- 
low Sn, except for the small NLDs of !!!Cd and !°'Dy at 
low excitation energy. The smaller result is due to the larger 
energy spacing of the theoretically calculated single-particle 
levels near the Fermi level. Overall, the results obtained by 
the CDFT combinatorial method are reliable. 


IV. SUMMARY AND PROSPECTS 


The combinatorial method has been adopted to describe the 
nuclear level densities for nuclear reaction calculations. The 
particle-hole state density is obtained with a combinatorial 
method using the single-particle level provided by the CDFT 
based PK1 effective interaction. After accounting of col- 
lective effects, including vibration and rotation, the energy-, 
spin-, and parity-dependent NLDs are obtained. Our results 
are compared with those from other NLD models, including 
phenomenological, microstatistical as well as non-relativistic 
HFB combinatorial models. The comparison suggests that 
besides some small deviation from different NLD models, the 
general trends among these models are basically the same. In 
conclusion, the CDFT combination method is capable to re- 
produce the experimental data at or below the neutron separa- 
tion energy. It implies that the CDFT combinatorial method 
is as reliable as other models to describe NLDs. Finally, the 
NLDs of the CDFT combinatorial method with normaliza- 
tion are compared with experimental data, obtaining excellent 
agreement with the observed cumulative number of levels at 


N(E) 


E (MeV) 


Fig. 7. (Color online) Comparison of the results obtained from CDFT combinatorial method (red lines) with the cumulative number of 
observed levels [49] (black lines) as a function of the excitation energy. 


373 low excitation energy and the measured NLDs below the neu- 
374 tron separation energy. 

35 Our results exhibit the predictive power of the CDFT com- 
binatorial method, even though the pairing correlations are 
not considered while the collective effects are empirical. In 
our following work, the CDFT combinatorial method can be 
improved by considering the inclusion of energy-dependent 
pairing correlations and taking a partition function approach 
381 to the treatment of vibrational enhancement. These results 
32 Will surely help our study on important neutron capture pro- 
cess such as r-process. 
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Fig. 8. (Color online) Comparison between the calculated NLDs based on the CDFT combinatorial method (red lines) and the experimental 
data [51, 52, 59-61, 64-68] (black dot). The full asterisk corresponds to the experimental data at the neutron separation energy Sn [49]. 
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